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1 Introduction 



Depending on the parameters involved, the dynamics of the film of viscous fluid can be described 
by different asymptotic equations. Under the assumption that the film is thin enough for viscous 
entrainment to compete with gravity, the time evolution model of a thin film of liquid on the inner 
surface of a cylinder rotating in a gravitational field was based on the lubrication approximation 
and examined by Benilov, O'Brien, and Sazonov [21 E]- The related Cauchy problem has the 
following form: 

yt + l[y]=0, y(0,x)=y o , y(-7t,t) = y(n,t), x G [ 7r, 7r] , t>0 (1.1) 

where 

d ( dy(x)\ 
l[y\ = — I (1 — a cosx)y(x) + b smx ■ — 1, a, b > (1.2) 

CL Ob V (X CC J 

Eigenmode solutions are very important in stability analysis, because even a single growing 
mode can destabilize an otherwise stable system. In case when all modes are bounded in time 
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and the corresponding eigenfunctions form a complete set, the system normally regarded as 
a stable one. Because, an arbitrary initial condition can be represented as a series of these 
eigenmodes; and since all of them are stable, so expected to be the solution to the initial- value 
problem. 

There are however counterexamples to the arguments above when each term of the series is 
bounded but the series as a whole diverges and the solution develops a singularity in a finite 
time. This effect was observed by Benilov, O'Brien, and Sazonov [2] for the problem 11.11 when 
parameter in (11.21) a = 0. For this case when the effect of gravitational drainage was neglected 
because of infinitesimally thin film they studied stability of the problem asymptotically and 
numerically. It was shown that even for infinitely smooth initial values numerical solutions blow 
up after a small number of iterations. 

The spectrum of the linear operator L that is defined by the operation and periodic 
boundary conditions y{—^) = y{^) for the special case when the parameter a = was stud- 
ied rigorously in [SI El E]- Using different approaches they justified that if the parameter b 
restricted to the interval [0, 2] then the operator L is well defined in the sense that it admits 
closure in L 2 (— 7r,7r) with non-empty resolvent set without breaking the boundary conditions 
y(—ii) = y{^)- The spectrum of the operator L is discrete and consists of simple pure imaginary 
eigenvalues only. As a result all eigenfunctions have the following symmetry y x (—x) = y\(x). 
The more general operator with the function sin(x) replaced by the arbitrary 27r-periodic func- 
tions was studied in pE] and it was proved that this operator multiplied by % belongs to a wide 
class of PT-symmetric operators which are not similar to self-adjoint but nevertheless possesses 
purely real spectrum due to some obvious and hidden symmetries. 

The phenomenon of the coexistence of the neutrally stable modes with explosive instability 
of the numerical solutions [2J (which correspond to drops of fluid forming on the ceiling of the 
cylinder where the effect of the gravity is the strongest) was studied analytically and explained 
in terms of the absence of the Riesz basis property of the set of eigenfunctions in [5]. The 
question of a conditional basis property of the set of eigenfunction is still open. 

For the case when a ^ 0, as it was discussed in [3J, the spectral properties of the operator L 
are not expected to differ a lot from the case a = 0. 

The goal of this paper is to find the factorization of the operator L (under some restrictions 
on parameters a and b) that would be in some sense similar to one we constructed for the special 
case a = in [7] (in this case the operator L is J-self- adjoint with the operator J defined as a shift 
J(f(x)) = f(jr — x)) and to examine some properties of the operator L using this factorization. 
The main difficulty to overcome here is an existing coupling between two subspaces spanned by 
positive and negative Fourier exponents which are not invariant subspaces of the operator L if 
a ^ 0. We also prove that the non- self-adjoint differential operator L has compact resolvent and 
as result spectrum of L is discrete with the only accumulating point at infinity. 
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2 Factorization of the non- self- adjoint operator L. 



We denote by D(T) and 9t(T) the domain and the range of linear operator T respectively. The 
notation C 2 is used for the standard Lebesgue space of scalar functions defined on the interval 
(— 7T, 7r). From here on L is the indefinite convection-diffusion operator L : 

'2 ^2 /r„.\/_\ rf \../-_\ , d 2/0*0 



£ I— ► £ , (Ly)(x) = —— [ (1 — a cosx)y(x) + b sinx 



<ix \ dx 

with the domain of all absolutely continuous 27r-periodic functions y(x) such that (Ly)(x) G £ 2 . 
In addition, we define the operator S : 

C 2 ^C 2 , (Sy)(x)=y'(x), 

where y'(x) G C 2 , y(— n) = y(n), and the operator M : C 2 i— > £ 2 , 

dy(x) 



(My)(x): = (1 — a cosx)y(x) + 5 sinx 



dx 



with the domain of all functions y(x) G £ 2 absolutely continuous on (—it, 0) U (0, 71") and such 
that (My)(x) G C 2 . Note, that, for example, y(x) = x -1 / 3 G D(M). The operator M can also 
be represented by the following expression 

(My) (x) = (1 — (a + b) cos x)y(x) + b (sin x • y(x))'. 

Theorem 1 If the parameters a and b satisfy the inequality 2a + b < 2, then L is a closed 
operator with a closed range and L = SM . 

Proof. Let us consider the operator A : 

C 2 ^C 2 , (Ay)(x) = (sin(x)y(x))' 
with D(A) = {y(x) \ y(x), (Ay)(x) G C 2 }. Then a function y(x) can be written as 

y{x) = — Lr • (c + f 0(t)dt), 0(t) G £ 2 . (2.1) 
sm(x) J 

If x > then 

0(t)dt| < ■ a(x) ■ x 1/2 , 

sin(x) 

where a(x) = ( \9(x)\ 2 dx) 1 ' 2 . Since the two summands in (12. ip have different orders of growth 
as x — > this implies that if y(x) G C 2 then c = and 



Moreover, 



|y(x)| (2.3) 
sm{x) 

A small modification of the same reasoning leads to the following estimation for every x G (— ir, n) 



2^|l/2 



with a(x) = | Jq \9(x)\ 2 dx\ 

Alternatively the same function y(x) can be written as 

y(x) = — ^ • (c - / 0(t)dt), G £ 2 . (2.5) 

with the same 0(x) as in ( 12. ip . Representation ( 12. 5ft yields the following relations 

y (a;) = ^L. f (2.6) 

and 



[y(x)|< 1 . f\ -P{x) (2.7) 



sin(x) 



with /3 (re) = |/;|^(x)| 2 rfx| 1/2 . 

It follows from and flZS]) that 



7T 



9{t)dt = 0. (2i 



o 



Starting from the point — it one can also obtain that 



y(x) = ■ 1 9(t)dt, (2.9) 



sin (x) 



\y(x)\ < { * { - 7 (x) (2.10) 
| sm(x)| 



with 7(x) = ( |6'(x)| 2 (ia:| 1/ ' 2 and 



o 



9{t)dt = 0. (2.11) 



Now we are ready to calculate M*. Using smooth functions y(x) such that y(x) = in some 
neighborhoods of the points — ir, and n (neighborhoods depend of y(x)) it easy to show that 

(M*z)(x) = (1 — a cosx)z(x) — b (sinx • z(x))' 
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for every z(x) G D(M*). Since the condition z(x) G T)(M*) yields 

z(x) G C 2 and (M*z){x) G C 2 , (2.12) 

for are fulfilled the conditions of the type (I2.4p . (12.71) and (12.101) . Taking into account the 
latter one can check that 

(My,z) = (y,M*z) 

for every y(x) G D(M) and z(x) under Conditions (12121) . Thus, D(M) = D(AT). The same 
reasoning shows that M " = M, so M is closed. 

Let 

(1 — a cosx)y(x) + b sinx • ^- - = m(x), y(x),u(x) G £ 2 . 
Our aim is to express y(x) via m(x). Let x G (— 7T, n), x ^ 0. Then ?/(a;) = 

(c(l + Signs)/2 + ci(l - Signa;)/2) ■ (sin |x|) a/fe ■ (cot |a;|/2) 1/b + 

^(sin|x|) a/fe - (cotlxl/2) 1 / 6 I* u(t)(sm\t\)-% ■ (sint)' 1 • (tan ltl/2) 1 ^, 



where c and ci are constants. The estimations that follow closely depend of a relation between 
a and b. We assume that 2a + 6 < 2. Then for x > 



M(i)(sint)- ( t +1 ) • (tant/2) 1/b dt = / v(t)(t)-^ +1) ■ (t/) 1/b dt, 

Jo 

where v(t) = M (t)(sint)-(t +1 ) • (tan*/2)V6(i)(f +1) • (t)" 1 / 6 , so, 

M (t)(sint)- ( t +1) • (tant/2) 1/b dt| < 



y=J=.^. (J r Wt )i^. 



Thus, the first summand (if c 7^ 0) for 7/(2;) has the order and the second one has the order 
x^ 1 ^ 2 a(x) with lima(x) = 0. Since u(;r) G £ 2 , c = 0. The same reasoning shows that c\ = 0. 

Thus, 

2/(ar) = ^(sin \x\) a/b ■ (cot |a;|/2) 1 / 6 / w(t)(sin |t|)~f • (sint)" 1 ■ (tan It^) 1 / 6 ^. (2.13) 

In particular, for m(x) e 1 we have 

y (x): = hsin \x\) a/b ■ (cot \x\/2) 1/b /"(sin ■ (sint)" 1 ■ (tan |t|/2) 1/6 di. 

</o 
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Some elementary estimations show that there are finite limits lim yo(x), lim yo(x) and lim yo(x) 

x^O x^—n x^tt 

with lim yo(x) = \imy (x). Let us show these relations. First, for t > we define w(t): = 

X — ► — 7T X — >7T 

(sst)-f- 1 . ^M±!2±fl\ Then Hm «;(*) = (1/2) 1 / 6 . Moreover, for x > 

y (x):= ^(sina;) a/fc • (cotrr/2) 1 / 6 f w{t)r^ 1+1 / b dt, 
b Jo 

so 

Ito(x) = i(sinx) a / fc • ( C0 tx/2f b -^w^ x )(xrf +1 / b , 
b 1 — a 

where £ x G (0, x). The latter yields y(0): = lim = y^-. Second, for t < ir we define w + (t): = 

(g)-f-i.((Z[z*). tan ( t /2)) 1/b . Then t lim w+(t) = 1 and for z(x): = (1 + a) -y Q (x) ■ ( sin(a;/2))* 
we have 

= i±^ (sina; )^ . f x w+a^-ty^-^dt. 

b Jo 

Let us fix e > Then there is 5 > such that 1 — e < w+(x) < 1 + e for every x G (it — 5, it). 
Next, for the same x 

z (x) = i^(sina;)^ • ( f $ w + (t)(>K - t)-^ 1 -^ b dt+ 
b Jo 

r w + (t)(u -t)- a *- i - i/b dt). 

Jn-5 



Since 
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with v XjS G (it - 5,7r), 



[ X w + (t)(n - ty^-^dt) = J^ w (v x , s )((ir - xy 1 ^ - 5-^) 
Jtt-s i + a 



b 



1 + a 



(I - e)((7r - x)- 1 ^ -5-^) < [ X w + (t)(n - tyi-^dt) < 



1 + a 

Moreover, for fixed 5 

!; . idJ / 

1 sill Tib ■ / 

t^ + TT-0 b 

and 



b /-, x// \ 1 + a t- i+£ \ 

(1 + e)((7r - x) — - 5 —). 



ii m l±^( sina; )^. r 5 w+ ( t )(n-ty^- 1 - 1 / b dt = o 



lim (sinx) & •tu(V 3 . < 5)o' * =0, 
so the equality \imy (x) = y^- is practically evident. Note also that yo(^) is even. 
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Thus, the function yo(x) is continuous on [— n, ir] and satisfies the periodic conditions. 

Now let u{x) = c + J* (j)(t)dt, where c = const and <f>(x) G C? . Then for x > 0: | L (p{t)dt\ < 
x 1 / 2 ( [n \4>(t)\ 2 ) . The latter estimation and (12.131) yield \imy(x) = c ■ yn(0). The same 
function can be re-written as following u(x) = c + + J ' <fi(t)dt or u(x) = c_ + J"* (p{t)dt. Then 
the estimations | j* <f){t)dt\ < (vr - x) x / 2 ( - /* |0(t)| 2 ) 1/2 , ar G (o, vr) and | j^(j){t)dt\ < (vr + 



,r 



|1 ^ 2 '' I^WI 2 ) 1 ^ 2 ) x ( —7T ) °) together with Representation (12. 131) yield lim y(x) = c + ■ ynfa 



X—*TT 



and lim y(x) = C--yn(—n). Thus, y(x) satisfies the periodic condition if and only u(x) satisfies. 

X— > — 7T 

The latter yields the equality 

L = SM. (2.14) 

Moreover, we shown that for every u(x) = c + f <f)(t)dt with (j)(x) G C 2 there is absolutely 
continuos on [— tt, tt] function y(x), such that (My)(x) = u(x), so $H(M) is dense in C 2 or, 
equivalently, Ker(M*) = {0}. 

Note, that M is boundedly invertible. Indeed, M = D + iC, where C: C 2 \— > C 2 , 

(Cy)(x): = z j- cosxy(x) — b (sinx ■ y{x))' 
and D: C 2 i-> £ 2 , 

{Dy){x)\ = (1 - (a+ -) cosx)?/(x) . 

It easy to check that C is self-adjoint. Moreover, .D is positive, bounded and boundedly invert- 
ible, so the problem of invertibility of M is equivalent to the problem of regularity of non-real 
numbers for a self-adjoint operator (for a more detail reasoning see, for instance, [7]). 

Now let us prove that L is closed. The operator S restricted on the subspace C\ C C 2 of 
functions orthogonal to constants has a bounded inverse. Let us find M~ l (Li). If (My)(x) G 
Ci, then .(My)(x)dx = 0, but f* n (y(x) smx)'dx = 0, so y(x) G M -1 (£i) if and only if 

y(x) G S(M) and (1 — (a + 6) cosx)y(x)c?x = 0. Let C 2 - = |(1 — (a + 6) cosx)| . Since 

for yo(x) we have 2n = (My )(x)dx = y (x)(l — (a + b) cosx)dx, yn(x) G" C 2 and 
C 2 = £ 2 +{jt/ ■ yo(x)}^c. Now let a sequence {y k (x)} be such that y k (x) — ► and ^(x) = 
(Ly k )(x) z(x). Then (My k )(x) = v k (x) + c fe , where w fc (x) = ((S , | £l )" 1 2:)(x), c fc = const, 
fc = 1,2, ... . Since (SI^J -1 is bounded, the sequence {v k (x)} has a limite v(x). In turn, in 
virtue of similar reasons the sequence w k (x) = (M~ 1 v k )(x) also has a limite. Simultaneously 
yk{x) = (M~ l {y k + C k ))(x) = w k (x) + c k yo(x). Thus, the sequence {c k } has a limite. The rest 
is straightforward. □ 

Corollary 1 If the parameters a and b satisfy the inequality 2a + b < 2, then the set D(L) is 
the linear sub-manifold H of the Sobolev space H l {— tt, n): 

®(L) = H:fe H\-7T, tt), /(tt) = /(-vr, ) sm(x)f G H\-n, vr) 

and is a Hilbert space with the norm defined as: 

12 - ll <%, + \\Mx)f'(x)\\ 2 Hl . 
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The reasoning of this corollary is the same as the reasoning of the corresponding proposition in 

For the next step we need the following simple remark. 

Lemma 1 Let H, be a Hilbert space, xi,X2 G 7i, (%i,Xz) 7^ 0. Let Tii. = {xi}- 1 , 7i 2 : = {x 2 } ■ 
Let P\ and P 2 be ortho-projections onto the subspaces Ji\ and 7i 2 respectively. Then P 2 \hi 
one-to-one mapping onto 7i 2 . 

Proof. Let 7i 3 : = Tii H 7i 2 . With no loss of generality we can assume that ||xi|| = ||x 2 || = 1, 
(xi,x 2 ) = a > 0. Then Hi — {fi- (x 2 - a ■ a;i)} MeC © H3, H 2 = {n ■ {x\ - a ■ x 2 )]>€C © ?4- Since 
= Ih 3 an d ^2(^2 — ot ■ %i) — a • (xx — a ■ x 2 ), the rest is evident. □ 

Theorem 2 If the parameters a and b satisfy the inequality 2a + b < 2, then the resolvent of L 
is compact of the Hilbert- Schmidt type. 



Proof. Let Co C C 2 be the subspace of constants, C±: = Cq. Since 
has the following matrix representation 



Ci, the operator L 









£11 



with respect to the decomposition C 2 = C © C±. The operator L w : 1 — > a ■ sinx is bounded, 
so we need to analyze the properties of L\\. From Theorem 1 we have L\\ = SM\c 1 . Let 
C 2 = M(D(M) n £1), y (x) = (M _1 l)(x), z (x) = ((M*)~ 1 l)(x). Then for y(x) E ©(M) n £1 
we have = (y, 1) = (y, M*z ) = (My,z ), so C 2 = {^o} -1 . From the other hand, (l,z ) = 
(My ,z ) = (y ,M*z ) = (yo,l). Since (see the proof of Theorem [T]) (y , 1) 7^ 0, the pair 
{l,^o} is under the conditions of Lemma [TJ Let Pi be the ortho-projection onto C\. Then 
£11 = S- {P x \c 2 ) ■ (M| £l ), so Ln 1 = (M" 1 !^) ■ (P1UJ- 1 • (S'UJ -1 - Thus ; L n is an operator of 
the Hilbert-Schmidt type. Since 



Rx{L) 



Rx(0) 
-R\(Lii)L w R\(L U ) 



the rest is straightforward. 
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